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Explicit Expressions for the Variance and Higher Moments of the Size of 
a Simnltaneous Core Partition and its Limiting Distribution 

By Shalosh B. EKHAD and Down ZEILBERGER 


Dedicated to William Y. C. ‘Bill’’ Chen, the tireless apostle of enumerative and algebraic com¬ 
binatorics in China (and beyond) 

Important Update (Sept. 1, 2015): It turns out that the second challenge below (except 
for the page limit) has been met before our paper was written, by Paul Johnson. See insightful 
comments by Marko Thiel and Nathan Williams: 

http://www.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/stcoreFeedback.html 

Note, in particular, that all our of theorems are now rigorously proved theorems. As mentioned 
by them. Theorems 2 and 3 have been anticipated in their paper “Strange Expectations” 

http://arxiv.org/abs/1502.07934 

A donation to the OEIS in honor of Paul Johnson, Marko Thiel, and Nathan Williams, has been 
made. The first challenge remains wide open! 

VERY IMPORTANT 

As in all our joint papers, the main point is not the article, but the accompanying Maple package, 
stCore, that may be downloaded, free of charge, from the webpage of this article 

http://www.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/stcore.html , 

where the readers can also find sample input and output files, that they are welcome to extend 
using their own computers. 

Introduction 

Many (perhaps most) combinatorial statistics (e.g. the number of Heads in tossing a coin n times, 
the number of inversions of an n-permutation [and, more generally, the number of occurrences 
of any pattern in an n-permutation]), are asymptotically normal, which means that if you denote 
it by Xn, figure out the average, a„ := E[Vji], and then figure out the variance, let’s call it 
m 2 (n) := E[(Xjj — a^)^], then the centralized and standardized version := {Xn — an)( 
tends, as n —oo, to the good old normal distribution (aka Caussian distribution) whose probability 
density function is -^=e~‘^ Our favorite way ([Z1][Z2]) of proving this is automatically, by using 
symbol-crunching to compute (at least the leading terms) of the general moment mr{n) and then 
prove, {automatically, of course) that lim„_>.oo equals 0 for r odd, and r!/((r/2)! 2”/^), for 

r even, the famous moments of the normal distribution. 
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But there are numerous exceptions! The most notable is the statistics “length of the longest 
increasing subsequence” defined over the set of n-permutations, where the intriguing Tracy-Widom 
distribution shows up. Another example is the subject matter of the present article, the random 
variable “size”, defined over the set of (s,t)-core partitions (see below for definition), where s and 
t are relatively prime positive integers. One of us (DZ) pledges a 100 dollars donation to the OEIS 
Foundation, in honor of the first prover(s), for meeting the following challenge. 

First Challenge: Prove, rigorously, that the scaled limiting distribution (see below) (as (s,t) both 
go to infinity, and s — t is a hxed constant) of the combinatorial random variable “size” defined on 
the set of (s,t)-core partitions is given by the continuous random variable 


OO 


E 




where Zk and Zk are jointly independent sequences of independent standard normal random vari¬ 
ables. 

One of us (SBF) verified that the first nine (standardized) moments of “size of an (s,t)-core parti¬ 
tion”, (as (s, t) —7- oo) converge to the corresponding moments of the above continuous distribution, 
and it is virtually certain that this is true in general (see below for details). 

This distribution is mentioned in [DGP], eq. (2.4), where it is called ?7ys(l). 

But the main purpose of this article is to show the power of symbol-crunching in deriving new 
mathematical knowledge. We will state deep new (polynomial) expressions for the first six moments 
of the random variable size of an (s, t)-core partition, and the first nine moments for the special case 
of (s, s -\- l)-core partitions. From the “religious-fanatical” viewpoint of the current “mainstream” 
mathematician, they are “just” conjectures, but nevertheless, they are absolutely certain (well, 
at least as absolutely certain as most proved theorems). We also briefly indicate how we derived 
these expressions, and indicate, for those obtuse mathematicians who would like to see mathematical 
proofs, how they may possibly be proved. 

We are not offering any prizes for just any old proofs, but we will be delighted to donate another 
$100 to the OFIS Foundation for meeting the following challenge. 

Second Challenge: Come up with a “hand-waving” and “soft” (yet rigorous!) a priori reason 
(whose length is not to exceed two pages) why the average, and any finite moment, must be a 
polynomial in s and t, and also come up with a “soft” (but rigorous!) upper bound for the degrees. 

This would, in one stroke, prove, rigorously, all the theorems in this article, because it would 
rigorously justify the (empirically) obvious fact that they belong to the polynomial ansatz, and 
hence discoverable by undetermined coefficients, by gathering enough data, and solving a (usually 
large) system of linear equations, thereby turning the undetermined coefficients of the desired 
polynomiai expressions into determined ones. (See below for details.) 
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Note that in many cases in combinatorics (e.g. [BZ], [E], [Z3]), such arguments are very easy, but 
in the present case, we don’t see it. We hope that one of our readers will! 

(s,t)-Core Partitions and Drew Armstrong’s Ex-Conjecture 

Recall that a partition is a non-increasing sequence of positive integers A = (Ai,..., Xk) with /c > 0, 
called its number of parts; n := Ai A^ is called its size, and we say that A is a partition ofn. 

Also recall that the Ferrers diagram (or equivalently, using empty squares rather than dots, Young 
diagram) of a partition A is obtained by placing, in a left-justified way, Xi dots at the i-th row. For 
example, the Ferrers diagram of the partition (5,4, 2,1,1) is 

* * * * 

* * 

* 


Recall also that the hook length of a dot {i,j) in the Ferrers diagram, 1 < j < A^, is the number of 
dots to its right (in the same row) plus the number of dots below it (in the same column) plus one 
(for itself), in other words A—i-|-A'— where X' is the conjugate partition, obtained by reversing 
the roles of rows and columns. (For example if A = (5,4, 2,1,1) as above, then A' = (5,3, 2, 2,1)). 

Here is a table of hook-lengths of the above partition, (5,4, 2,1,1): 

9 6 4 3 1 

7 4 11 
4 1 

2 
1 

It follows that its set of hook-lengths is {1, 2,3,4,6, 7,9}. A partition is called an s-core if none of 
its hook-lengths is s. For example, the above partition, (5,4, 2,1,1), is a 5-core, and an i-core for 
all i > 10. 

A partition is a simultaneous {s, t)-coie partition if it avoids both s and t. For example the above 
partition, (5,4,2,1,1), is a (5, ll)-core partition (and a (5,12)-core partition, and a (100,103)-core 
partition etc.). 

For a lucid and engaging account, see [AHJ]. 

As mentioned in [AHJ], Jaclyn Anderson ([A]) very elegantly proved the following. 

Theorem 0: If s and t are relatively prime positive integers, then there are exactly 

{s -\-t- 1)! 
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(s, t)-core partitions. 

For example, here are the (3 + 5 — l)!/(3!5!) = 7 (3,5)-core partitions: 

{empty, 1,2,11,31,211,4211} . 


Drew Armstrong ([AHJ], conjecture 2.6) conjectured, what is now the following theorem. 

Theorem 1: The average size of an (s,t)-core partition is given by the nice polynomial 

(s — l)(t — l)(s + t + 1) 

24 


For example, the (respective) sizes of the above-mentioned (3, 5)-core partitions are 

0,1,2,2,4,4,8 , 

hence the average size is 

0 + l + 2 + 2 + 4: + 4 + 8 21 ^ 

7 ~ Y ~ 

and this agrees with Armstrong’s conjecture, since 

(3-l)(5-l)(3 + 5 + l) 


24 


= 3 


Armstrong’s conjecture was recently proved by Paul Johnson ([J]) using a very complicated (but 
ingenious!) argument (that does much more). Shortly after, and almost simultaneously (no pun 
intended) it was re-proved by Victor Wang [Wan], using another ingenious (and even more com¬ 
plicated) argument, that also does much more, in particular, proving an intriguing conjecture of 
Tewodros Amdeberhan and Emily Sergei Leven ([AL]). Prior to the full proofs by Johnson and 
Wang, Richard Stanley and Fabrizio Zanello [SZ] came up with a nice (but rather ad hoc) proof of 
the important special case of (s, s + l)-core partitions. 

We should also mention the very interesting approach in [YZZ], that proved an important special 
case (also proved in [AL]) of a conjecture of Amdeberhan and Leven (but not directly related to 
the subject matter of the present article). 

Higher Moments 

Recall that the r-th moment about the mean of a random variable, X, is E[(V — E[V])’'], where 
E is the expectation operation. 

Theorem 2: If s and t are relatively prime positive integers, then the variance (aka as the second 
moment (about the mean)) of the random variable “size of an (s,t)-core partition”, is given by the 
nice polynomial expression: 

--^st{t-l){s-l){s + t + l){s + t) . 
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For example for (3, 5)-core partitions it equals 


^3-5(5-1) (3-1) (3 + 5 + 1) (3 + 5) = 6 
and indeed (since the average is 3) 


(0 - 3)2 + (1 - 3)2 + (2 - 3)2 + (2 - 3)2 + (4 - 3)2 + (4 - 3)2 + (8 - 3)2 


Theorem 3: If s and t are relatively prime positive integers, then the third moment (about the 
mean) of the random variable “size of an (s,t)-core partition” is given by the nice polynomial 
expression: 

st {t — 1) {s — 1) {s + t + 1) (s + t) (2 s‘^t + 2 st"^ — 3 s‘^ — 3 st — 3t‘^ — 3) 

For example for (3, 5)-core partitions it equals 

1 90 

3-5-(5 - l)-(3 - l)-(3 + 5 + l)-(3 + 5)-(2 • 32 • 5 + 2 • 3 • 52 - 3 • 32 - 3 • 3 • 5 - 3 • 5^ - 3) = — 
60480 7 

and indeed (since the average is 3) 

(0 - 3)3 + (1 - 3)3 + (2 - 3)3 + (2 - 3)3 + (4 - 3)3 + (4 - 3)3 + (8 - 3)3 _ 90 

7 “ y 


Theorem 4: If s and t are relatively prime positive integers, then the fourth moment (about the 
mean) of the random variable “size of an (s,t)-core partition” is given by the nice polynomial 
expression: 


1 


st {t — 1) {s — 1) {s + t + 1) {s + t) 


4838400 

(19 + 38 s^t^ + 19 s'^t'^ - 51 sH - 102 s^t‘^ - 102 s'^t^ - 51 st'^ + 36 

+72 s^t + 108 s2t2 72 st^ + 36t‘^ - 33 ^2^ - 33 st"^ + 36 ^2 + 36 st + 36t^ + 120) . 


For example for (3,5)-core partitions it equals and indeed (since the average is 3) 


(0 - 3)4 + (1 - 3)4 + (2 - 3)4 + (2 - 3)4 + (4 - 3)4 + (4 - 3)4 + (8 - 3)4 726 

7 “ ~ 
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Theorem 5; If s and t are relatively prime positive integers, then the fifth moment (about the 
mean) of the random variable “size of an (s,t)-core partition” is given by the nice polynomial 
expression (in computerish): 

l/95800320*s*t*(t-l)*(s-l)*(s+t+l)*(s+t)*(46*s**6*t**3+138*s**5*t**4+138*s**4*t**5+46 
*s*=t'3*t=t'=t'6 -211*s**6*t**2 - 633*s*=t'5*t=t'*3~844*s*=t'4*t=t=*4 - 633*s=t'=t'3*t**5-211*s**2*t**6 
+ 333*s**6*t + 999*s**5*t**2 + 1665*s**4*t**3+1665*s**3*t**4 + 999*s**2*t**5+333*s*t**6 

- 180*s**6-540*s**5*t - 1283*s**4*t**2-1666*s**3*t**3 - 1283*s**2*t**4 - 540*s*t**5 
-180*t**6 + 420*s**4*t + 840*s**3*t**2 + 840*s**2*t**3 + 420*s*t**4-180*s**4 - 360*s**3*t 

- 540*s**2*t**2 - 360*s*t**3 - 180*t**4 + 327*s**2*t + 327*s*t**2 - 180*s**2 - 180*s*t 

- 180*t**2 - 3780) 

For example for (3,5)-core partitions it equals and indeed (since the average is 3) 

(0 - 3)5 + (1 - 3)5 + (2 - 3)5 + (2 - 3)5 + (4 - 3)5 + (4 - 3)5 + (8 - 3)5 2850 

7 “ ~Y~ 

Theorem 6: If s and t are relatively prime positive integers, then the sixth moment (about the 
mean) of the random variable “size of an (s,t)-core partition” is given by the nice polynomial 
expression (in computerish); 

l/4184557977600*s*t*(t-1)*(s-1)*(s+t+1)*(s+t)*(307561*s**8*t**4+1230244*s**7*t**5+ 
1845366*s**6*t**6+1230244*s**5*t**7+30756l*s**4*t**8-2056306*3**8*t**3-8225224*s**7*t**4- 
14394142*s**6*t**5-14394142*s**5*t**6-8225224*3**4*t**7-2056306*s**3*t**8+537206l*s**8*t 
**2+21488244*3**7*t**3+42976488*3**6*t**4+53720610*3**5*t**5+42976488*3**4*t**6+21488244 
*s**3*t**7+5372061*3**2*t**8-6453396*3**8*t-25813584*3**7*t**2-60704054*3**6*t**3- 
91764618*3**5*t**4-91764618*3**4*t**5-60704054*3**3*t**6-25813584*3**2*t**7-6453396*3*t 
**8+2985120*3**8+11940480*3**7*t+39743142*3**6*t**2+77437746*3**5*t**3+96285048*3**4*t** 
4+77437746*3**3*t**5+39743142*3**2*t**6+11940480*3*t**7+2985120*t**8-11104272*3**6*t- 
33312816*3**5*t**2-55521360*3**4*t**3-55521360*3**3*t**4-33312816*3**2*t**5-l1104272*3* 
t**6+2985120*3**6+8955360*3**5*t+23840061*3**4*t**2+32754522*3**3*t**3+23840061*3**2*t** 
4+8955360*3*t**5+2985120*t**6-9109476*3**4*t-18218952*3**3*t**2-18218952*3**2*t**3- 
9109476*3*t**4+2985120*3**4+5970240*3**3*t+8955360*3**2*t**2+5970240*3*t**3+2985120*t 
**4+8664840*3**2*t+8664840*3*t**2-62687520*3**2-62687520*3*t-62687520*t**2+626875200) 


For example for (3, 5)-core partitions it equals 2346, and indeed (since the average is 3) 

(0 - 3)5 + (1 - 3)5 + (2 - 3)5 + (2 - 3)5 + (4 - 3)5 + (4 - 3)® + (8 - 3)® 16422 

---= —-— = 2346 . 

7 7 

The last three theorems regard the special case of (s, s + l)-core partitions. 

Theorem 7; If s is a positive integer, then the seventh moment (about the mean) of the ran¬ 
dom variable “size of an (s, s -|- l)-core partition” is given by the nice polynomial expression (in 
computerish): 
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l/149448499200*s**2*(s-l)*(s-2)*(2*s+l)* (124496*s**14 - 527660*s**13 - 127268*s**12+ 
2133077*5**11 + 1565655*s**10-3928575*s**9-7848989*s**8-3573289*s**7 + 7257797*5**6 
+16741975*5**5+16528197*5**4+3583272*5**3 - 67819248*5**2-18541440*5+138620160) 

*(5+1)**2 

Theorem 8: If s is a positive integer, then the 8-t/i moment (about the mean) of the random variable 
“size of an (s, s + l)-core partition” is given by the nice polynomial expression (in computerish): 

1/914624815104000*5**2 * ( 5 -I)*(2*5+1)* (308851624*5**18 - 2759073420*5**17 + 7345195650*5**16 
+ 1614779679*5**15 - 27716691813*5**14 - 3203324556*5**13 + 61922226136*5**12 + 52270343442*5*^ 

- 49025878614*5**10 - 146716496688*5**9 - 153171599682*5**8 - 30342055161*5**7 + 158893451131*: 

- 165853921776*5**5 + 1073038790016*5**4 + 9260929255680*5**3 - 11293714925568*5**2 

- 19188060088320*5 + 21924617379840)* (5+l)**2 

Theorem 9: If s is a positive integer, then the ninth moment (about the mean) of the random 
variable “size of an (s, s + l)-core partition” is given by the nice polynomial expression (in com¬ 
puterish): 

1/182467650613248000*5**2*(5-1)*(5-2)*(2*5+1)* (28092743584*5**20 - 284614603048*5**19 
+ 908242721124*5**18 - 87722680542*5**17 - 4040707469643*5**16 + 1347179583168*5**15 
+ 11350317109273*5**14 + 4824122583716*5**13 - 15816684214230*5**12 - 31535118689736*5**11 

- 29475404073738*5**10 + 2671156715274*5**9 + 63014451511513*5**8 + 79700408583680*5**7 

+ 45859575725901*5**6 - 377516262865248*5**5 + 6309067352294376*5**4 + 10737857697068736*5**3 

- 38301852570773760*5**2 - 26103018295756800*5 + 48704747653094400 )* (5+l)**2 

A Crash course in Combinatorial Statistics 

Recall ([Z1][Z2]) that given a sequence of combinatorial random variables (e.g. the number of 
Heads upon tossing a fair coin n times), whose combinatorial generating function is either known 
explicitly (C'n(t) = {l + t)^ in this trivial case), or only as expressions in t, but for many n, one first 
finds the discrete probability generating funetion Pn(t) := Cn{t)/Cn{l) (e.g. Pn{t) = (1 -|- t)‘^(2‘^ 
for coin-tossing), (under the uniform distribution). Of course H„(l) = 1 as it should. To find the 
first moment, avn (alias average, alias mean, alias expectation) one either derives explicitly H))(l) 

(n/2 in this trivial case), or “guesses” in more complicated cases. Next one centralizes getting 
the centralized probability generating function Pn{t) = P^(t)/t“”", and gets the variance (alias 
second moment about the mean), m 2 (n), by computing (or “guessing”) {t-^YPn{t)\t=i, and more 
generally, the higher moments, mk{n), are given by (t^)^Pn{t)\t=i- From these one derives the 
standardized (scaled) moments ak{n) := ■ It is often the case, as n —)• 00 , that the ak{n) 

tend to fixed numbers, in which case we have a limiting distribution. 

The above is easily extended to the case where the sequence of combinatorial random variables de¬ 
pends on several discrete parameters, like the present case where they depend on the two parameters 
s and t. 
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The Wonder of Internet Searches: How we found the (so far conjectured, but abso¬ 
lutely certain) Scaled Limiting Distribution 

As mentioned above, more often than not, this limiting distribution is the good-old Gaussian, but 
this is definitely not the case this time. Using our Theorems 2-4, we computed that the limiting 
skewness (standardized third moment) happens to be 

4 ;— 

lim 03 ( 5 , t) = - ■ VlO PS 1.807 , 

s,t—^oo 7 

while the limiting kurtosis (standardized fourth moment) turns out to be 

57 

lim 04 ( 5 , t) = -— ~ 8.1429 . 

s,t—>-oo 7 

Being citizens of our time, we googled 
1.807 8.1429 statistics 
and, lo and behold, got the link 

WWW.aueb.gr/conferences/Crete2015/Papers/Dalla.pdf , 

(reference [DGP]), that mentioned that these are the skewness and kurtosis of the continuous 
probability distribution that they call Uvs{^) (already mentioned in the introduction), but we will 
call Z 

°° ^2 I -2 

y ._ ~r 

■ ^ 47r2A:2 ’ 

fc=i 

where Zk and Zk are jointly independent sequences of independent normal random variables. 

Lindas Giraitis, one of the authors of [DGP], kindly offered the following interesting information 
via email. 


“VS” stands for “rescaled variance statistic”. Characterization of the U 1/5(1) distribution as the sum of the weighted 
sum of iid normals was well- known in statistical literature. It was established by Geoffrey S. Watson ([Wat]) in the 
context of goodness-of-fit tests on a circle, see formula ( 15 ) in his paper. Later in [GKLT], the authors show that 
the VS statistic has the same limit as the limit derived by Watson in the context of goodness-of-fit tests on a circle. 

A Human Intermezzo: A Crash Course in Moment Generating Functions 

This inspired us to compute higher moments of the continuous probability distribution, Z, as 
follows. 

Recall that the moment generating function of a probability distribution X is the (exponential) 
generating function 

00 

k=0 



It is well-known and trivial to see that, for any fixed constant, c, 

Mcx{t) = Mx{ct) , 

and it is also well-known and easy (but not utterly trivial) to see that if X and Y are independent 
random variables, then 

Mx+vit) = Mx{t)MY{t) . 

Hence if is a sequence of pairwise independent random variables and is a sequence 

of positive numbers such that their sum converges, then 

OO 

i=l 


If, in addition, the X^’s are identically distributed, denoting Mxiit) by M{t), we have 


i—\ 


The famous moments of the standard normal distribution, z, are, as mentioned above, 0 for r odd 
and r!/(2’'/^(r/2)!) for r even. Hence the r-th moment of is (2r)!/(2’’r!). Hence the (exponential) 
moment generating function is 


M^2(t) 


^(2r)! r y-(2r)!,, 

Y 2^t*! v\ ' 

r=0 * * r=0 


(1 - 2 t )“^/2 


Hence the (exponential) moment generating function of + 5^, where z and 5 are independent 
standard normal distributions, is 

M,2+J2(t) = ((1 - 2t)-'/2^' = (1 - 2t)-i . 

Hence the (exponential) moment generating function of the continuous distribution Z is 

OO 

fc=i 

But thanks to good old Leonhard’s iconic sin(7ra;) = vrxn^i(l “ f^-), this equals 

sin ^/t/2 

and this can be used to find as many (straight) moments as desired (take the coefficient of t^ 
in the Maclaurin expansion and multiply by kl). From these straight moments one easily com¬ 
putes the moments about the mean, using, thanks to the binomial theorem, E[(X — mi)^] = 
■*mi, and from them the standardized moments. 
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[This is implemented in procedure VSmoms(N) in the Maple package stCore]. 


Using this, one of us (SBE) found: 

4 


57 


ag = - . VTO , 04 = ^ 


820 /— 

05 = — ■ vTo , 


06 — 


1537805 

7007 


466860 /— 193032265 

«7 = • VlO , 08 = 


Og — 


1001 
70231858960 
2263261 


17017 


• Vw 


and these coincide exactly (up to the ninth moment) with the limiting (scaled) moments of our 
combinatorial random variable “size of an (s, t)-core partition”, as {s, t) go to infinity, (and s — t is 
bounded). 

How did we derive the above Explicit Expressions for the First Nine Moments? 

Let us now briefly describe how we were able to discover the above very deep theorems regarding 
the hrst nine moments. 


The hrst step was to use Anderson’s bijection, as modihed in [AHJ], and to work with (s,t)-Dyck 
Paths, that are two-dimensional walks, using horizontal and vertical positive unit steps, from the 
origin to the point {s,t) staying above the line sy — tx = 0. 

[Equivalently, the number of ways of tossing a coin t S times such that you win S dollars if it lands Heads and 
lose t dollars if it is Tails, breaking even at the end, but never being in debt]. 

Each such path has associated with it a certain set of positive integers, nicely described in [AHJ], 
and also important is how many such positive integers there are. Introducing the formal variables 
q and w, we associate the weight to each such path, where s is the sum of the labels and k is 
the number of such positive labels (see [AHJ] for exactly how to determine them). 

The same dehnitions apply to partial walks from (0,0) to (i,j) where we look at the vertical steps 
{i,j — 1) ^ (b j); and we naturally associate the partial weight due to each such vertical step 
(assuming {s,t) are fixed) as follows (for the sake of typographical clarity we use x * *y for x^). 

Wt{i,j){q, tc) := g * * sj — ti' — b \ - tc** 1 

; s j — t i' — b>0} J ; s j — t i'— b>0} 

then the weight-enumerator of the set of partial walks from (0,0) to (f,j) may be computed by the 
recurrence (we suppress the implied dependence on {q,w)) 

FsAhi) = Fs,t{i+Wt{i,j-1) ■ Fs,t{i,j-l) . 



We impose the obvious initial condition Fs^t{0,0) = 1, and boundary conditions Fs^t{i,j) = 0 if 
i < 0 or sj — ti < 0. 
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Using this linear recurrence scheme, one can compute these weight-enumerators very fast. 

But Fg^tii, is but a stepping stone. At the end of the day we take Fs^t{s,t){q,w), expand it 

in terms of powers of w and perform the umbral substitution 

yjk ^-k{k-l)/2 _ 

(This corresponds to going from the vector of hook-lengths of the first column of the (s, t)-core 
partition to its shape, see [AHJ]). 

All this is implemented in procedure Fab(a,b,q) in our Maple package stCore. 

Using this efficient, Dynamical Programming, approach we can crank owt these weight-enumerators, 
Fs^t{(l)i for many choices of coprime pairs {s,t), and by repeatedly applying the operator q-^, 
and then plugging-in q = 1 collect numerical data for the respective moments, and from these, 
numerical data (for many distinct {s,t)) about the moments about the mean. Then, inspired by 
Drew Armstrong’s ex-conjecture, that expressed the average as a polynomial, we assume (without 
any a priori theoretical justification!, see the second challenge) that these are polynomials (i.e. we 
use the polynomial ansatz), and use undetermined coefficients, and fit the data with polynomials. 

A possible Approach for (most probably ugly) mathematical proofs 

The difficulty with the Fg^tiQ), in their initial definition as weight-enumerators, according to “size 
of (s,t)-core partitions”, is that they are only defined for s and t relatively prime. While the 
notion of (s, t)-Dyck path makes sense when s and t are not relatively prime, the nice formula that 
enumerates them, (s 1 — l)!/(s!t!), is no longer valid. 

A natural approach would be not to be hung up on “nice” formulas, and set-up a general recursion 
schemes for the weight-enumerators of partial walks, Fs^t{i,j){(l,w), but no longer restricted to 
relatively prime s, t. Then get some recurrence scheme that would enable an (ugly, ad-hoc) proof 
in the so-called holonomic ansatz. But since this approach is only likely to produce ugly proofs, 
and we are absolutely certain that our Theorems 2-9 are correct, we rather not even try, and do 
not encourage anyone to follow this approach. (But those who insist should start with the easier 
special case when t = s + 1.) 

But we will be thrilled if auy of our readers would meet our two challeuges stated iu 
the iutroductiou! 
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